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f-{— ' ! Abstract 

c/3 ■ In this paper we analyze the effect of a combined pure homogeneous strain 



and simple shear in a principal plane of the latter on the propagation of sur- 
face waves for an incompressible isotropic elastic half-space whose boundary 
is normal to the glide planes of the shear. This generalizes previous work 

'^ . in which, separately, pure homogeneous strain and simple shear were con- 

sidered. For a special class of materials the secular equation is obtained in 
O ' explicit form and then specialized to recover results obtained previously for 

the two cases mentioned above. A method for obtaining the secular equa- 
tion for a general form of strain-energy function is then outlined. In general 

J^ I this is very lengthy and the result is not listed, but, for the case in which 

Q«^ ' there is no normal stress on the half-space boundary, the result is given, for 

illustration, in respect of the so-called generalized Varga material. Numer- 
ical results are given to show how the surface wave speed depends on both 

^^ ■ the underlying pure homogeneous strain and the superimposed simple shear. 

*yN . Further numerical results are provided for the Gent model of limiting chain 

extensibility. 

>: 
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1 Introduction 

The vibration and wave propagation characteristics of rubberhke sohds are im- 
portant for practical apphcations in, for example, engine mountings and seismic 
isolators, and an understanding of these characteristics is therefore of crucial im- 
portance. There have been many contributions in recent years to the analysis of, in 
particular, vibrations and surface and interfacial waves in finitely deformed elastic 
bodies, but relatively little attention has been given to shearing deformations, for 
which the orientation of the principal axes of deformation depends on the magni- 
tude of the shear. In this paper we are concerned with the influence of a simple 
shear combined with a pure homogeneous strain on the propagation of surface 
waves. 

Analysis of surface waves in a half-space subject to pure homogeneous defor- 
mation dates back to the classic paper of Hayes and Rivlin P, which was con- 
cerned primarily with compressible materials, while Flavin [2] considered the same 
problem for an incompressible material, specifically the Mooney-Rivlin material 
(for completeness, mention should be made of early and noteworthy articles by 
Biot |3] and by Buckens [5] and the related work of Biot summarized in [5J). A 
more detailed analysis was provided for general incompressible elastic materials by 
Dowaikh and Ogden [6], who also examined the connection between surface wave 
propagation and stability of the half-space, again for the case of pure homogeneous 
strain. We refer to this latter paper, that by Chadwick [7] and the recent contribu- 
tion by Fu |8j, for pointers to the relevant literature. The only paper thus far that 
deals with the influence of shear on surface wave propagation appears to be that 
by Connor and Ogden [9] , and the closely related paper [10] . These are concerned 
with simple shear and give explicit results for a certain class of incompressible 
materials, of which the best known representative is the Mooney-Rivlin model for 
rubberlike solids. 

The purpose of this paper is to extend the analysis in the above papers to the 
situation in which a simple shear is combined with a pure homogeneous deformation 
and for a general incompressible isotropic hyperelastic material. These extensions 
are motivated by several factors. First, one can easily picture situations of practical 
interest where a solid is stretched and sheared: for instance, a rubber isolator under 
a bridge is subjected to a vertical compression and is then sheared as a result of 
thermal extensions and contractions of the roadway. Also, the combination of 
stretch and shear theoretically allows a square block of incompressible material to 
be maintained in a deformed state with no shear tractions nor normal tractions 
on its faces [TT], in contrast to simple shear and its associated Poynting effect. 
Finally, the quite recent application of nonlinear elasticity to the modelling of some 
biological tissues (see Humphrey [12] or Holzapfel and Ogden [13]) has highlighted 



the need for solutions valid for more general strain-energy functions because the 
functions used for biomaterials are usually more complicated than those used for 
rubberlike materials. 

In Section 2 the basic equations for the (static) finite deformation consisting of 
a pure homogeneous strain (or triaxial stretch) followed by a simple shear in a prin- 
cipal plane of the pure homogeneous strain are summarized for an incompressible 
isotropic elastic material. The corresponding equations for infinitesimal surface 
waves superimposed on this finite deformation are then outlined. The equations 
and zero incremental traction boundary conditions are applied to surface waves 
propagating in the direction of shear and parallel to a half-space whose normal lies 
within the plane of shear. The secular equation for the wave speed is then obtained 
in implicit form in respect of a general strain-energy function. 

This equation in examined in detail in Section 3. For a class of materials that 
includes the Mooney-Rivlin strain energy, the secular equation is given in simple 
explicit form, which relates the wave speed to the material parameters, the de- 
formation and the normal stress on the half-space associated with the underlying 
finite deformation. Next, an alternative approach is outlined, which enables the 
secular equation to be made explicit in the more general case. However, since 
this is very complicated we illustrate the results by considering a generalized form 
of the Varga material model. Specifically, in the absence of normal loads on the 
half-space boundary we obtain the secular equation in relatively simple explicit 
form as a quartic for the squared wave speed. On the other hand, in the presence 
of normal loads we give the corresponding explicit form of bifurcation condition 
(corresponding to zero wave speed), which consists of two separate equations. For 
each of the material models considered numerical results are described in Section 
4, illustrating the dependence of the wave speed separately on the stretching and 
shearing parts of the underlying deformation. Finally, for comparison, some lim- 
ited results are given for a model of finite chain extensibility of rubberlike solids 
originally due to Gent [T^, and recently used by Horgan and Saccomandi [I5] to 
model strain-hardening arterial walls. 

In closing this Introduction, we note that the method and the results presented 
in the paper are easily extended to a generic type of homogeneous pre-strain having 
a plane of symmetry, and to unconstrained (compressible) solids. The assumptions 
of combined stretch and shear and of incompressibility lead to some simplifications 
and to neat algebraic expressions, but they are inessential and should not be viewed 
as restrictive. 
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(a) Unit cube (b) Triaxial stretch (c) Simple shear 

Figure 1: Triaxial stretch followed by simple shear of a unit cube. 

2 Basic equations 

2.1 Finite static deformations 

Let (Xi,X2,X3) be a rectangular Cartesian coordinate system. Consider a homo- 
geneous semi-infinite body occupying the region X2 > in its natural (unstressed) 
configuration. Suppose that the body is composed of an incompressible isotropic 
elastic material characterized by a constant mass density p and a strain-energy 
function W{\i, A2, A3) (per unit volume) written as a (symmetric) function of the 
(positive) principal stretches (Ai,A2,A3), which satisfy the incompressibility con- 
straint 

A1A2A3 = 1. (2.1) 

The body is subjected first to a finite pure homogeneous triaxial stretch with 
principal stretch ratios /ii,/i2,/i3 (/Wi/i2/i3 = 1) and then to a finite simple shear 
of amount k, so that a particle initially at (Xi,X2,X3) is displaced to (i:i,a;2,a;3) 
according to 

Xi = fiiXi + KP2X2, X2 = 112X2, £3 = /i3-^3- (2.2) 

Figure 1 shows the successive deformations of a unit cube in the half-space; in 
Figure 1(c) the angle is the angle of shear, which is given by tan0 = k. 

The deformation gradient F associated with the deformation 02. 2|) has compo- 
nents 

/il K.p2 
/i2 . (2.3) 

/i3 



Let Xi,X2,a;3, be the rectangular Cartesian coordinates associated with the (or- 
thogonal) principal axes of FF (the left Cauchy-Green deformation tensor), one 
of which is the £3 direction, so that X3 = £3. Within the (xi,X2) plane let the 
Xi principal axis make an angle 6 with the xi axis measured in the anticlock- 
wise sense. Then the coordinates {xi,X2,X3) are related to {xi,X2,Xs) through a 
rotation according to 



Xi 



^'ijXj'f [^'^jj 



COS 6 — sin 6 

sin 9 cos 9 

1 



(2.4) 



The principal stretches, subject to (12. ip . and the angle 9 are given in terms of 
the components (12. 3 p by 



Af cos' 9 + ^ sin' 9 

{\l — A2) cos 6^ sin 6^ = fx^K 

A^sin^e + A^cos^e 



2 I 22 



f^l, 



(2.5) 



from which we obtain 



Ai ± A2 = Y (/Ui ± 122^ + K,'^f4, A3 = l/(/ii/i2), 



(2.6) 



and 



tan 29 = 2/i^fi;/(/i? - f4 + K^f4) (2.7) 

(see, for example, [HI [HI [E]). When the material is sheared but not stretched 
(/xi = /i2 = /^3 = 1), for example, the well-known formulas 



Ai^2 = (V4 + fi;2±/t)/2, A3 = l, tan2e = 2/ 



K, 



(2i 



are recovered. 

Let (Ti, (72, (73 denote the principal Cauchy stresses. Since the material is isotropic 
the principal axes of Cauchy stress coincide with those of FF , and hence, analo- 
gously to (12.51) . we may express the components of the Cauchy stress on the axes 
associated with the coordinates {xi,X2,X3), denoted by ctjj, as 



o"ii = cTi cos^ 9 + a2 sin^ 9, 
<5"i2 = (o"i — (^2) sin 9 cos 9, 
o"22 = o"i sin^ ^ + o"2 cos^ 9, 



(2.9) 



with (Ti3 = 0-23 = 0, 0-33 = a3. 



In terms of the strain-energy function, the principal Cauchy stresses are given 

by 

dW 

aj = Aj— p, zG {1,2,3} no sum over z, (2.10) 

oXi 

where p is a Lagrange muhipher associated with the constrain (12. ip . 

The initial half-space X2 > becomes the half-space £2 > after the deforma- 
tion described above, and the load on its boundary £2 = has normal component 
(T22 and shear component 0"i2. We now consider the propagation of surface waves 
in this half-space. 

2.2 Superimposed infinitesimal surface waves 

Now consider a surface (Rayleigh) wave propagating with speed v and wave num- 
ber k in the Xi direction in the considered deformed half-space (£2 > 0), with 
attenuation in the £2 direction. We emphasize, as noted above and in contrast 
to many previous studies of surface waves in deformed materials, these directions 
do not coincide with the principal directions (xi,X2) of strain for the underlying 
pre-deformation. It seems that the only paper considering the propagation of non- 
principal surface waves over a non-principal plane is that by Connor and Ogden 
[9], where the problem was solved for a sheared half-space {fii = fi2 = fis = I) and 
for a specific class of materials (see also the related paper by Connor and Ogden 

uni). 

Let the incremental displacement u associated with the wave have components 
Ui, i G {1,2,3} in the coordinates system (a;i,X2,£3), with dependence, in gen- 
eral, on (xi,a;2,a;3, t), where t is time. The corresponding components Sij of the 
incremental nominal stress tensor are then given by 

^ dui dui 
Sij = Bijki-^^ + PgTr - pSij, (2.11) 

where Bijki are the components of B, the tensor of first-order elastic moduli as- 
sociated with the underlying deformation, and p is the incremental counterpart of 
p. The non- vanishing components of the elasticity tensor B, denoted Bij^i, in the 
coordinate system {xi,X2,X3) are given by 

Biijj = XiXjWij, 

S.,,, = (XiW, - X,W,)X^/iX^ - X'^) I ^ J, 

Bijji = Bjiij = Bijij — XiWi i ^ i (2-12) 

(see, for example, [El El [I9]), where Wi := dW/dXi, Wij := d'^W/dXidXj. Note 
that expressions for Bijij in the case Xi = Xj, i ^ j, may be obtained by a limiting 
process. 



The two sets of components of B are related by 

and they satisfy the symmetry properties 

Bijki = Bkiij, Bijki = Bkiij- (2-14) 

Also, note that for i, j G {1, 2}, in particular, B has six non-zero entries in the 
(xi,X2,X3) system, namely Bun, ^2222, ^1122, ^1221, ^1212, and^ ^2121,^ whereas 
in the {xi,X2,X'i) system it has ten: 5iiii,^52222,^Sii22,^ -B1221, 5m2, ^2121, the 
counterparts of those above, together with -B1112, -B1121, -B1222, and -82122- 

For surface waves propagating on a half-space subject to a pure homogeneous 
stretch and no shear (k = 0), Dowaikh and Ogden [H] found it convenient to 
introduce the shorthand notations 

« = -B1212, 1 = B2121 = \i \a, 2/3 = i?iiii + i?2222 — 2-Bii22 — 2-Bi22i- (2.15) 

The corresponding quantities in the (xi, £2, £3) system are obtained by using fl2.13p 
in the form 



5i2i2 = a cos^ e + 2P sin^ 9 cos^ ^ + 7 sin^ 9, 
B2121 = a sin^ 9 + 2/3 sin^ 9 cos^ ^ + 7 cos^ 9, 
Bun + ^2222 - 2^1122 - 2£i22i = 2/3 - 6(2/3 - a - 7) sin^ 9 cos" 9. (2.16) 



It is easily checked that neither a, 7, /3 nor d, 7, /3 depend separately on the 
constant Cauchy stress components (cti, (72, Cs) (or equivalently on aij). It follows 
from (EUD that 

d-7 = (a-7)cos20, d + 7 - 2/3 = (a + 7 - 2/3) cos4e. (2.17) 

From the latter it follows that if the equality a + 7 — 2/3 = is satisfied, as it is for 
certain forms of the strain-energy function, then we also have d + 7 — 2/^ = for 
any angle 9. In such a case the equality is said to be 'structurally invariant' [2U] . 
Finally in this section, we note that the strong ellipticity condition may be 
expressed in the form 

a > 0, /3 + yd7>0 (2.18) 

(see, for example, Dowaikh and Ogden [6] and Chadwick [7J). It follows immedi- 
ately that d > 0, but the counterpart of the second inequality in ( I2.18p for the 
{xi,X2,X3) system is not so simple, and is not therefore given here. 



2.3 Equations of motion and boundary conditions 

We now restrict attention to waves of the form 

u(xi,X2,t) = U(e2)e^'=(^^-''*\ (2.19) 

independent of £3, where U is a function of kx2, k is the wave number and v 
the wave speed. This must satisfy the incremental equations of motion and the 
incremental incompressibility constraint. Respectively, these are 

^ = P^, u,. = ina;2>0. (2.20) 

Additionally, we impose the condition of zero incremental traction on the plane 
boundary, so that 

S2» = onx2 = 0, 2 = 1,2,3. (2.21) 

In view of equations (I2.19p and (12. lip . Sij and p take forms similar to that of 
the displacement u, that is 

{%(xi,X2,t),p(xi,X2,t)} = {kS,,{kx2):ikP{kx2)V^^'''-'"\ (2.22) 

where we have introduced the functions Sij and P of kx2 It is a straightforward 
exercise to check that the in-plane equations decouple from the out-of-plane equa- 
tions, so that there is no loss of generality in taking f/3 = 0. The incremental 
equations of motion and incompressibility (l2.2Up and (I2.2ip then specialize to 

i5n + 5^i = -Xf/i, iSi2 + S'^^ = -XU2, if/i + t/^ = 0, (2.23) 

where X := pv^ and a prime signifies differentiation with respect to the argument 

kX2. 

It is convenient to introduce the quantities 

n = -i5'2i, T2 = -15*22, (2.24) 

which are functions of £2 proportional to the incremental tractions on the planes 
parallel to the boundary. Thus, from the incremental boundary conditions (I2.2ip . 
we have 

ri(0) = r2(0) = 0. (2.25) 

The expressions 

Su = KBnn + pWi + i5ni2f/2 + Bu2iU[ + Bn22U^ - iP, 

ir2 = ii?ii22f/i + i-Bi222f^2 + B2i22U[ + (-B2222 + P)^2 ^ i-^' 

in = iBiuiUi + i{Bu2i + p)U2 + lU[ + 52122?^^, 

S12 = iBnuUi + iaU2 + (^1221 + p)U[ + Bu22U^ (2.26) 

8 



may be derived from equations f l2.19p . fl2.22l) and (12.111) . 
The short-hand notations defined by 



hi 



-S1222 
-B1121 



-B1112 
-B2122 



I [7 
I [7 



a 
a 



{213 ~a--f) cos 29] sm2e, 
{2(3 -a--f)cos2e]sm2e, 



(2.27) 



prove useful in rewriting the equations of motion in a convenient form. Indeed, 
a few manipulations and use of fl2.23p . fl2.26|) and f l2.27p lead to the first-order 
differential system 

^' = iNt |:=[f/i,f/2,ri,r2]^, 



(2.28) 



where 
N 



Ni Na' 



Ni 



-z/21/7 -(7-'3-22)/7 
-1 



N, 



1/7 0' 




K(i) 



X - 2(/3 + 7 - 0-22) + z>|i/7 z>i2 + i>2i(7 - '3-22)/7 
1^12 + '>2i(7- 0-22 )/7 X - d + (7 - 0-22)77 



(2.29) 



The classical approach to solving the equations of motion (12.281) subject to 
(I2.25P consists of writing the solution as a decaying exponential function, such as 

^{kx2) = ^^e"i^''\ %q) > 0, (2.30) 

and then solving the eigenvalue problem N^ = g^ . Its solution is of the form 

^{kx2) = Ail^e'^^^^^ + v42^2e'^2'=^^ (2.31) 

where Ai, A2 are constants such that (12.251) is satisfied. Explicitly, gi and g2 are 
the complex roots of the propagation condition det(N — ql) = with positive 
imaginary parts, where I is the fourth-order identity. The resulting equation is the 
quartic 

7g' + 2z>2ig^ + (2/3 - X)q^ + 2uuq + a - X = (2.32) 

for q, while ^\ i = 1,2 may be taken as proportional to any column vector of the 
matrix adjoint to N — gjl. If we take the third such column, for example, then 

C = [Qi, -Qi, QiilQi + i'2iqi - 7 + (3-22), -(7 - '3-22)g,^ + z>i2gi + a- X]^ . (2.33) 

The boundary conditions ^(0) = [[/i(0), ^^2(0), 0, 0]"^ then yield a homogeneous 
system of two linear equations for the two unknowns Ai, A2, whose determinant 
must be zero. This leads to the secular equation for the surface wave speed, written 
in the implicit form 



7(« - ^){<li + (I1Q2 + 92) - (7 - ^22){iqiq2 +a-X) 

+ ['>i2'>2i - (7 - o-22)^]gig2 + [vi2iqiq2 + «>2i(tt - X)]{qi + ^2) 
It remains implicit as long as gi and g2 are not expressed explicitly. 



0. (2.34) 



3 Explicit secular equations 

Here the problem of an incremental surface wave propagating over the surface 
of a stretched and sheared incompressible half-space is solved, up to the explicit 
derivation of the secular equation giving the speed of the Rayleigh wave. 

3.1 Secular equation for materials with 2/? = a + 7 

Consider the Mooney-Rivlin strain-energy function 

W = \C{\\ + Ai + A2 - 3) + \D{XlXl + XlXl + A^A? - 3), (3.1) 

where C and D are constant material parameters. For (13. ip . a, 7 and /3 reduce to 

a = Xl{C + DXl), ^ = Xl{C + DXl), 2^ = {X\ + Xl){C + DXl), (3.2) 

and so it is apparent that Mooney-Rivlin materials belong to the general class of 
materials with energy function such that 

2/3 = a + 7. (3.3) 

This class has the remarkable property that the propagation condition for inhomo- 
geneous plane waves with displacement proportional to 3fJ{exp[i/c(a;i -|- qx2 — vt)]} 
always admits g^ + 1 as a factor when the body is maintained in an arbitrary state 
of static homogeneous deformation. In fact, Pichugin [21J proved that this class is 
the most general one admitting this factorization. 

Returning to the propagation condition (12.321) and using (12.51) . we find that 
when (13. 3p holds the quantities appearing in the quartic (12.321) simplify to 

" = iK'^ (/^i + lAi<^) > 7 = 7A^VL ^^21 = i>i2 = -iX^'^iiIk, 

2P = ^Xf[^Jil + ^^l{l + ^^)], X = ^Xffil{fil/fil-r]'), (3.4) 

where the quantity 

V .= ^Jl-4- (3.5) 

/^2 V /^i 7 

has been introduced in the last equality. Then, after removal of the factor 7 7^ 0, 
the quartic (12.321) factorizes as expected to give 

{q^ + l){q^ - 2Kq + K^ + 7]^) = 0, (3.6) 

and its roots with positive imaginary parts are 

gi = i, g2 = « + ^V- (3-7) 
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In passing, we note that the condition 3(gj) > for decay of the wave yields an 
upper bound on the admissible wave speed, such that 

0<X = pv^ < fijX^h, (3.8) 

with the limiting value /i^A^^7 corresponding to a body wave not a surface wave 
(which requires strict inequality). 

Finally, substituting the values (13. 7p for gi and g2 into (I2.34p . we arrive at the 
explicit secular equation 

7]"^ + 7]^ + {3 + K^ - 2a22)v - (1 - ^22)' = 0, (3.9) 

where CT22 := ^2^22 / (fJ^'il) ■ This secular equation is identical (apart from notational 
differences) to that obtained when the material is stretched but not sheared (k = 0) 
[21 E] or when the material is sheared but not stretched {pi = 1) [S]. 

In the case of the Mooney-Rivlin strain energy (13. ip . 7A^^ = C + D / {H1112Y , so 
that both r] and o"22 are independent of the amount of shear k. Of course rjR, the 
positive real root of the cubic (13. 9p . depends on k. By differentiation we find that 
driR/d^K^) = —Vr/I'^Vr + Vr + {^ ~ ^22)^] < and so that, by (13. 4p ^. the surface 
wave speed is a monotone increasing function of the amount of shear. Hence for 
Mooney-Rivlin materials subject to a given triaxial stretch /ii,/i2,/i3, the greater 
the subsequent shear is, the faster the Rayleigh wave propagates. 

In the limit X = the secular equation (13. 9p reduces to the bifurcation equation. 
This defines a curve in (/ij, n, a22)-space that separates the region where the surface 
wave exists (and is unique) and the half-space is stable from the region where 
the surface wave does not exist and the half-space might become unstable (see 
Chadwick and Jarvis [22j or Dowaikh and Ogden [6j). The bifurcation equation is 

^y + f ^ V + (3 + «:^ - 2^22) (^) - (1 - a22)^ = 0. (3.10) 

P2J V/^2/ V/^2/ 

When there is no normal Cauchy pre-stress in the £2 direction {0^22 = 0), the bifur- 
cation criterion becomes universal with respect to the entire class of incompressible 
materials satisfying (13. 3p since it does not depend on the material parameters as- 
sociated with a specific strain energy, but only on the pre-stretch and the amount 
of shear. 

3.2 Secular equation for a general strain energy 

For a general (incompressible, isotropic) strain-energy function (2/3 7^ a -|- 7), the 
propagation condition (I2.32p does not factorize as in (13. 6 p and it is difficult in 
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general to use the quartic (12.321) for the derivation of an exphcit secular equation 
Instead, we employ the method of the 'polarization vector' devised by Currie 
and based on the fundamental equations 



U(0) ■ [k(")u(o)] = 0. 



(3.11) 



Here, the symmetric 2x2 matrix K*^") is the lower left block of N", where n is 
any positive or negative integer. Thus, for instance, K*^^) is given by (I2.29p ^. and 
we find that K*^^^) and K*^^^ are given by 



X-a 



7(X - a) + (7 - 6-22)^ 
-z>2i(X -a) + z>i2(7 - 0-22) 



ly: 



12 



-Z>2l(X -a) + Z>12(7 - Cr22) 

(X - a)(X -2/3 -27 + 2(722 ) 



(3.12) 



and 



Ki) 



'(1) 



2i)2iK\l> /^ + 2K- 



'(1) 



(1) 
12 



Ki) 



'(1) 



Z>2li^i^77 + ^22 + (7 - '^22)i^liV7 



(1) 



i>2li^l277 + ^^2^ + (7 - 5-22)i^liV7 



(1) 



2(7 - ^22)i^l277 



(1) 



(2) (2) (2) 

respectively. Here Kl^^ and KI2 are linear in X and 7^22 is independent of X 

Now, the fundamental equations (13.111) . written in turn for n 
the homogeneous system 



(3.13) 
X. 
1,-1,2, yield 



-"■11 
-"■11 

'(2) 



Kll' K 



-"■12 
-'^12 

(2) 



-"■22 
-"■22 



12 



K. 



(2) 
22 



f/l(0)f/2(0)+C.C. 
f/2(0)f/2(0) 



(3.14) 



where c.c. denotes the complex conjugate of the preceding term. The determinant 
of the left-hand matrix in (I3.14p must be zero, and the resulting equation is the 
explicit secular equation for surface waves in a stretched and sheared incompressible 
material. It is a quartic in X, but since its explicit expression is very lengthy we 
do not give it here, although it has been obtained using Maple and Mathematica. 
We note, however, that in the stress-free configuration ((J22 = 0,/ij = 1, k = 0) it 
factorizes into the product of a term linear in X and the cubic of Rayleigh [23| for 
isotropic incompressible elastic solids, namely 



C=^-8C' + 24C-16 = 0, C--=pv''/l. 



(3.15) 



which has a unique real root C,r ~ 0.9126. Here 7 > is the shear modulus of the 
incompressible material in that configuration. 

As an illustration, we consider the generalized Varga strain-energy function 



W = C(Ai + A2 + A3 - 3) + D(AiA2 + A2A3 + A3A1 - 3), 



(3.16) 
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where again C and D are constant material parameters, but different from those 
in (13. II) . For discussions of fl3.16p see, for example, Carroll [25], Haughton [2S], or 
Hill [27]. For fl3A6|) . a, 7 and /3 reduce to 

^_ XliC + DXs) ^_ XliC + DXs) ^ _ X^X2iC + DXs) ^^ ^^^ 

Xi + A2 Xi -\- X2 Xi + A2 

and it is then apparent that f l3.16p belongs to the class of materials with energy 
function such that 

/3 = ^07. (3.18) 

For such materials we find first that, on use of (12.51) . 

" = X2(x +A )2 ^^^ ^ ^^^^ ^ /^2'^^)^ 7 = pp;— ^yy (/^i/^2 + /i2)^ 
Aal^Ai + A2) 

^^12 = -2 ,2/x , X N2 /^2^(^l + /^l/^2 + f^W), 
A2[Ai + A2J 

^^21 = -2-2TT— — --^/i2«(/^l/^2 + /i2)- (3-19) 

A2[Al + A2J 

Secondly, we find that, in the absence of normal load ((T22 = 0) on the plane 
boundary in the underlying configuration, the secular equation simplifies to 

2(/ii + fi2)x'^ - (3/ii + 7;U2)/i^a;^ + [(/ii + 5/i2)/i^ + 2fi2{3fil + 5fil)]^^x'^ 

-4(/ii + /i2)^/i2/i^°a; + (3/ii - /i2)(/ii + /U2)^Ai2/^^^ = 0, (3.20) 

where 

x:= A^(Ai + A2)^X/7, /i2:=(/ii + /i2)2 + /i^/€l (3.21) 

The bifurcation equation becomes simply 

3;Ui - /i2 = 0, (3.22) 

which is universal for the class of materials satisfying (13.171) . and is independent 
of the amount of shear, i.e. it is the same bifurcation equation as for a material 
subject only to a triaxial stretch [2E]- When the material is stretched but not 
sheared (k = 0), the secular equation factorizes into the product of a term linear 
in X and the cubic in the secular equation for principal surface waves |.6j. When 
the material is sheared but not stretched (/ij = 1), on the other hand, it reduces 
to 

2a;^ - 5/iV + (8 + 3/i2)/a;2 - 8/°a; + V^ = 0, /x^ = 4 + K^ (3.23) 
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with X appropriately specialized. 

Finally, we find that, in the presence of normal load (o"22 7^ 0) on the plane 
boundary in the underlying configuration, the secular equation is too lengthy to 
reproduce here, but the bifurcation equation is relatively simple and consists of the 
two possibilities 

a22 + /i2(3/il -yU2)[(/il +/i2)^ + /U2'«^] = 0> 0^22 = /iV2(/Ul + 7^2) , (3.24) 

where 0^22 = A2(Ai + X'lfo^ij^- For stability 0^22 must lie between the values 
specified by these two equations. This result generalizes that in Dowaikh and 
Ogden [S] (their equation (7.15)), which is applicable for k = 0, and recovers their 
result in this case. 



4 Numerical results 

4.1 Mooney-Rivlin materials 

For materials such as the neo-Hookean or Mooney-Rivlin models satisfying 2/^ = 
a + 7, the secular equation (13.91) for stretched and sheared bodies has been written 
in a form which is formally identical to that obtained by Connor and Ogden [10] 
for sheared (but not stretched) bodies, with the difference that their r] and 0^22 
are a specialization of ours to the case /ii = /i2 = /^s = 1- As a result the 
ijli'^ri) dispersion curves are common to those in [TU] and there is little interest in 
reproducing them here. 

Instead we concentrate on the bifurcation condition (13.101) . For a Mooney- 
Rivlin material subject to a pure homogeneous deformation (k = 0) and no normal 
load ((T22 = 0), Biot [5] found that (linearized) surface instability occurs at a critical 
stretch of (/Ui)cr = 0.544 for the plane strain /i3 = 1 and of (/ii)cr = 0.444 for the 
equibiaxial pre-strain fj,2 = fJ^s, and Green and Zerna ^29j found a critical stretch 
of (/ii)cr = 0.666 for the equibiaxial pre-strain /xi = fi^. For a stretched (/Xj 7^ 1) 
and sheared (k 7^ 0) Mooney-Rivlin material subject to a pre-load (J22 7^ we find 
from (13.101) that /ii//i2 is a monotonic decreasing function of the amount of shear 
K. Hence, in the cases where /ii//i2 is a monotonic increasing function of /ii, we 
conclude that if the half-space is stable in the stretched configuration, then it will 
remain stable in the stretched and sheared configuration whatever the amounts of 
pre-load and of shear are applied. These cases include the plane strain ^3 = 1 
(because then fii/fi2 = fA) ^^^ the equibiaxial pre-strains ^2 = A^s, A^i = A^s 
(because then /ii//i2 = /i/ ,/ii, respectively). Figures 2(a) and 2(b) show indeed 
this behaviour for (/ii)cr in the plane strain case, for different values of the scaled 
normal load 0^22. 
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Figure 2: Plot of (/ii)cr as a function of k, for a Mooney-Rivlin material in 
plane strain under the normal loads (from top curve to bottom curve) (a) 0^22 = 
1.5, 1.75, 2, 2.25, 2.5, 2.75, 3; and (b) 0^22 = -2.5, -2, -1.5, -1, -0.5, (thick curve) 
0.5. 



4.2 Varga materials 

Now we turn our attention to a half-space made of generalized Varga material 
03.16p . subject to large stretches and shear such that no extension occurs in the 
direction normal to the plane of shear. Hence /is = A3 = 1 and /i2 = yu]~^. 
When there is no normal load (o"22 = 0), then the critical stretch found from 
(I3.22P is (/Xi)cr = 0.577 [28j. As /ii is increased we find that the wave speed 
first increases and then decreases, whatever the amount of superposed shear, with 
Cs := \/{C + D)/{2p), the speed of shear waves in the undeformed material, as an 
asymptotic limit. In contrast to Mooney-Rivlin materials, the presence of shear 
decreases the surface wave speed in Varga materials subject to a given (fixed) 
amount of triaxial stretch /ii,/i2,/i3. Figure 3(a) illustrates these features, with a 
plot of the squared surface wave speed, scaled with respect to cj., and traced as a 
function of /ii for different values of k. At k = 0, /ii = 1, the material is unstressed 
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and we recover f^/c^ = 0.9126, as indicated by the dotted lines. We record here 
that for some values of the parameters equation (13.201) may have two solutions for 
X. However, because of the rationalization implicit in the derivation of fl3.20p . the 
second solution is spurious and is not associated with a surface wave. 

Next we exploit the remarkable property of the deformation (12. 2 p that it allows 
a square block made of a hyperelastic material to be stretched and sheared without 
either shear tractions or normal tractions applied on its faces [TT]. This possibility 
arises when fif = {1 + k^)ii\. The actual dependence of the /Xj on n is then found 
from the constitutive equation of a given material. For the neo-Hookean model, 
Rajagopal and Wineman [Hj found that 

/ii = (l + fi:2)^, /,2 = /is = (1 + /t^)"i (4.1) 

It can be checked that the same condition applies for Mooney-Rivlin and for gen- 
erahzed Varga materials. We computed the speed of a surface wave travelling over 
such a block, made of a classical Varga material (D = in (I3.16P ). Figure 3(b) 
displays the variation of the squared wave speed, scaled with respect to the squared 
shear wave speed c^ = C/(2p) in the undeformed material, as a function of the 
amount of shear. As k increases, the /ij vary according to (14.11) . and the wave speed 
increases in a monotone manner, at least within the 'reasonable' range [0, 5] for n. 
For greater, but unrealistic, values of n we found that the wave speed v eventually 
decreases, with c^ as an asymptotic limit. 

4.3 Sheared Gent materials 

We conclude with the study of the strain-energy function proposed by Gent [H], 
namely 

W = \CJr^\n[l-^^±^^±^^^^, A? + A^ + A1<3 + J^, (4.2) 

that accounts for the limiting chain extensibility of some elastic materials such as 
rubber or soft biological tissue. Here C is the infinitesimal shear modulus and Jm 
is a constant. Recently, Horgan and Saccomandi [15] showed how and to what 
extent this model could be used to describe finite deformations of strain-stiffening 
biological tissues, such as aortas and arteries. Typically, these biomaterials can 
be subjected to quite large deformations (within certain limitations) and exhibit 
highly nonlinear constitutive equations. These features are quite well accounted 
for by the strain-energy function (14. 2p : it is nonlinear in the stretches, the limiting 
condition (I4.2p 9 imposes an upper bound on the value of the stretch, and the 
two parameters C and J^ can be adjusted to render a satisfactory picture of 
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Figure 3: Plot of f^/Cg (a) against fii for (from top curve to bottom curve) k = 
0, 1, 2, 3, and (b) against k in plane strain for a square block of Varga material. 

many deformations and of the strain-hardening process (see [15] and the references 
therein to other works by Horgan and collaborators for extensive discussion of and 
justification for the use of the Gent strain energy-function for the modelling of 
stress-hardening biological tissues). Using classical experimental data ^U\ on the 
aorta of a 21 year old male and on the aorta of a 70 year old male, Horgan and 
Saccomandi computed values of Jm = 2.289, 0.422 for the younger and older aorta, 
respectively. 

Here we focus on a body composed of Gent material subject to a finite shear 
(yUj = 1, K > in (12. 2p ). We see at once from (2.7) that the limiting chain 
extensibility condition (I4.2p 9 yields an upper bond for the amount of shear and 
thus for the angle of shear. 



K 



\ J n 



< tan 



(4.3) 



Hence, using the values of Jm given above, we find that the younger aorta cannot 
be sheared beyond an amount of shear of approximately 1.513, corresponding to 
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an angle of shear of about 56.54°. For the older (stiffer) aorta, the corresponding 
values are: 0.650 and 33.01°, respectively. 

Now turning back to small amplitude motions superposed on a large homoge- 
neous deformation with principal stretches Ai, A2, A3, we find that for Gent mate- 
rials the elastic moduli (2.14) are given by 

^ — n ™ \2 ^, _ ^ "^m ,2 

t« — ^ ~i ~ To To To'^l) / — ^~^ ~ To To Tn^T 



2/3 = C 



<^m + 3 — A]^ — A2 



\2' 



Jm + 3 — A^ — A2 



X2 



A? + A^ + 



<^m + 3 — A^ — A2 — A3 

2(A^ - A^)^ 

J^Yn \ i — Aj^ — A2 — Aq 



(4.4) 



When the homogeneous deformation is a shear with amount k, they reduce to 



O O ri 



a. 



4 7 



K^ 



:(V4+^+fi:) 



G Jr, 



7 



4 J^ - fi;2 



(V4 



K^ — fi;)2, 



2/3 = C 



fJrr 



7 — K-2 
•J rn. f^ 



. 2/t2(4 + ^2 
2 + fi:2+ ^ 



Orr 



K^ 



(4.5) 



These expressions, together with the relations (12. 5p . evaluated for ^i = fi2 = 1, 
allow for the explicit computation of the quantities a, 7, /3 appearing in the secular 
equation. By (I2.16p . we have 



a 



C- 



Jr» 



J rr 



K^ 



k' + 



2k'' 



J rr 



K^ 



7 = C 



Jyi 



J rn. ^ 



1 + 



2k' 



J rn, '^ 



2/3 = C 



1^12, ^^21 



Uyy 



J ri 



K^ 



2 + k' 



2k^{2 - K^) 



-C 



Jfy 



J rn, '^ 



It 



•Jm ^ 

2k' 



(4.6) 



Now yj pv' /C, the surface wave speed scaled with respect to the speed of a (bulk) 
shear wave propagating in the undeformed material, can be plotted as a function 
of the amount of shear k G [O, \/ Jm ) • Figure 3 shows that the surface wave 
speed increases monotonically with the amount of shear, mildly for small shears 
and then dramatically as the limiting amount of shear is approached (vertical 
asymptotes at k = 1.513, 0.650 for the younger and older aorta, respectively). The 
curves are plotted in the absence of normal load (o"22 = 0); for a discussion on the 
influence of normal load on the wave speed and on the surface stability of sheared 
incompressible materials, see Connor and Ogden [9]. 
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